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Abstract 

In this work we analyze the effect of smoothing maps containing arrival directions of cosmic rays with a 
gaussian kernel and kernels of the mexican hat wavelet family of order 1, 2 and 3. The analysis is performed by 
calculating the amplification of the signal-to-noise ratio for several anisotropy patterns (noise) and different 
number of events coming from a simulated source (signal) for an ideal detector capable of observing the full 
sky with equal probability. We extend this analysis for a virtual detector located within the array of detectors 
of the Pierre Auger Observatory, considering an acceptance law. 



1 Introduction 

The origin, chemical composition and mechanisms of acceleration of the ultra-high energy cosmic rays (UHECRs) 
still remain a mistery half a century after their discovery [1], and are subject of great interest by the scientific 
community due to their intrinsic relation with particle physics, cosmology and astrophysics. The flux of cosmic 
rays with energies above 1 EeV (10 18 eV) is approximately 1 particle per square kilometer per year and 1 particle 
per square kilometer per century for energies above 60 EeV[5]. 

To achieve a reasonable statistic of events with energies above 1 EeV, it was built in the province of Mendoza, 
Argentina, the Pierre Auger Observatory, the result of an international effort of 17 countries. This observatory is 
a pioneer in the application of the hybrid technique of simultaneous fluorescence and surface detection. The first 
consists on the measurement of the fluorescence light resulting from the de-excitation of the nitrogen molecules 
which interacted with particles from the extensive air shower (EAS) generated by the primary particle. The 
surface detection consists on the detection of particles from the EAS by an array of 1,600 water-Cherenkov tanks 
composing a triangular grid with 1.5 km spacing from each other, covering an area of approximately 3,000 km 2 . 



The identification of possible astrophysical sources and the investigation of the magnetic fields which per- 
meates the universe are studied by analyzing the arrival directions of cosmic rays. The correlation of these 
directions with the large scale distribution of matter in the universe, such as the galactic and supergalactic 
planes are considered large scale anisotropics. A small scale anisotropy is characterized by the association be- 
tween arrival directions of cosmic rays and point sources, such as stars, distant galaxies and other kind of objects 
which are far enough to be considered point sources. 

Since the sources of UHECRs are still unknown, it is important to search for them. Amongst the candidates 
to sources of UHECRs are the active galactic nuclei, gamma-ray bursts and magnetars[3 . Here we analyze the 
performance of several convolution kernels on searching for anisotropies possibly correlating to some astrophysical 
sources of UHECRs. 



2 Wavelets 

Wavelets are mathematical functions belonging to the L 2 space that satisfy some requirements. The first one is 
the admissibility condition, 

— : — : — duo < oo, (1) 
M 

where tp(t) is the wavelet. This condition guarantees the reconstruction of a signal without loss of information. 
Also, the wavelet shall satisfy the condition of zero norm, i. e., 

i/>(t)dt = 0. (2) 

This last condition means that the average value of the wavelet in time domain shall be zero at zero frequency [3]. 
The continuous wavelet transform (CWT) may be formally written as: 

3(a,r)= [ f(tW T (t)dt, (3) 



where s (s > 0, s £ R) and r (r € R) are, respectively, the scale (dilation) and translation parameters. So, the 
CWT decomposes a function f(t) in a basis of wavelet $> StT (t). The inverse transform is given by: 

/(f) = J J $(s,T)V StT (t)dTds. (4) 

These wavelets ~$! aT (t) are obtained from a so-called mother-wavelet, by changing the dilation and translation 
parameters as follows: 

*.,r(«) = 4* ( — ) ■ (5) 
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The mexican hat wavelet family (MHWF) and its extension on the sphere have been widely used aiming the 
detection of point sources in maps of cosmic microwave background (CMB)[51 [7], due to the amplification of 
the signal-to-noise ratio when going from the real space to wavelet space. 

The MHWF is obtained by successive applications of the laplacian operator to the two-dimensional gaussian. 
A generic member of this family is: 

C— Tl™ 

= i^r v ( 6 ) 

where <j> is the two-dimensional gaussian (4>{x) = ^e~ x / 2<T ) and the laplacian operator is applied n times. 



3 Celestial Maps 

When studying anisotropies, a powerful tool are the celestial maps, which are segmented pixelizations of the 
celestial sphere. The events map is the celestial map which represents, in two dimensions, the arrival directions 
of cosmic ray events in the celestial sphere. 

Due to limitations of the detector itself, it is impossible to determine the exact arrival direction of an event. 
Each event detected is convolved with a probability distribution related to the angular resolution of the detector, 
that is, for each event there is an associated point spreading function (PSF). Therefore, it is extremely useful 
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to convolve the celestial maps with functions associated to the PSF of the detector, aiming to maximize the 
signal-to-noise ratio. Mathematically, this process of convolution (filteringM) may be written as 



M f (fb) =aj M(r)$(r, fh)dSl, (7) 

where a is a normalization constant, M(f) is the number of cosmic ray events in the direction r, $(r, r~o) is 
the kerneQ of the transformation and fo is the position vector representing each point in which the integral is 
evaluated. In the discrete case this process is: 

Mf[k)= Ej fJ) • (8) 

where M(J) is the number of cosmic rays associated to the pixel of index j in the direction fj. 



4 Analysis Procedure 

In the present work we have tested the capability of the MHWF to detect point sources embedded in different 
backgrounds. The simulated source was simulated in the position (I, &)=(320°, 30°) (galactic coordinates) in the 
sky, with a dispersion of 2°. We have considered three possible scenarios for the source, varying the amplitude 
of the source with respect to the background, and the number of events simulated in its direction. In the first 
case the source is 10% stronger than the background and 100, 000 events were simulated. In the second case, 
we have simulated 5, 000 events and the source has an amplitude of 1% with respect to the background. In the 
third case there are only 50 events coming from the direction of the source, which has an amplitude of 0.1% 
with respect to the background. 

For the background we have simulated 500, 000 events according to several anisotropy patterns. These 
anisotropy patterns were divided in two classes: with and without the acceptance of the detector. In the case 
without acceptance, the distribution of events in the sky is uniform, modulated only by the anisotropy pattern 
imposed by the simulation. In the case with acceptance, we considered a detector in a place with latitude 35° 
28' 00" S and longitude 65° 18' 41" W (approximately the coordinates of the Pierre Auger Observatory). We 
have assumed that the event are detected according to a zenith angle distribution that follows sinOcosO, where 
0° < 9 < 60° is the zenith angle[S]. 

Four anisotropy patterns were simulated as the background maps, with and without the acceptance of the 
detector (figures [l] to [4}. In all the cases these 500, 000 events were simulated. The simulated anisotropy patterns 
were: 

• isotropic: isotropic distribution of events; 

• dipole 1: a dipole with excess in the galactic center (l,b) = (0°,0°), with amplitude 7% with respect to 
the background; 

• dipole 2: a dipole with excess in the direction (l,b) = (266.5°, —29°), with amplitude 0.5% with respect 
to the background; 

• sources: several sources with different angular scales a and amplitudes A in the directions (l,b): (0°,0°) 
[a = 7.0°, A = 100%], (320°, 90°) [a = 1.5°, A = 5%], (320°, -40°) [a = 0.5°, A = 1%], (220°, 10°) 
[a = 3.0°, A = 5%], (100°, -70°) [a = 2°, A = 10%], (240°, 50°) [a = 20°, A = 5%], (350°, -80°) 
[a = 6.0°, A = 0.5%], (100°, 50°) [a = 30°, A = 50%], (140°, -40°) [a = 4.0°, A = 200%] and (60°, 50°) 
[a = 3.0°, A = 2%]. 

In order to verify the power of identification of the wavelets, we calculated the amplification (A) of the 
signal-to-noise ratio, which is: 

A=^, (9) 

where wq is the value of the central pixel associated to the source in the non-filtered source map, uif is the value 
of the same pixel in the filtered source map, <7o is the root mean square (RMS) of the non-filtered background 
map and <jf is the RMS of the filtered background map. 



1 In this paper we make no distinction between the processes of convolution, smoothing and filtering. 
2 We also make no distinction between the terms kernel and filter. 
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Figure 3: Simulated background maps with (left) and without (right) acceptance, for the dipole 2 case. 



According to Gonzalez-Nuevo et al [5], the maximum amplification of the signal-to-noise ratio for a source 
with dispersion 70 embedded in a white nois^J background is obtained by convolving the maps with a gaussian 
kernel with dispersion 7 = 70- This shall be used as a reference to check the consistency of our results. 

3 White noise has the property of being homogeneous and isotropic, with a uniform power spectrum. 
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Simulated Events Map containing several sources (with acceptance law) 



Simulated Events Map containing several sources (full sky) 




Figure 4: Simulated background maps with (left) and without (right) acceptance, for the sources case. 



5 Results and Discussion 

Using equation [9] we have calculated the amplifications in all the simulated cases, with scales spanning from 0.2° 
up to 6.0°. We present here only the graphs of amplification in which the source has an amplitude of 1% with 
respect to the background, with 5,000 events in its direction. The results are shown in figures [5] to [8j 




Figure 5: Amplification for the isotropic case with (left) and without (right) acceptance. 

We are interested in obtaining amplifications greater than one, indicating that the signal-to-noise ratio from 
the filtered map is greater than the one from the non-filtered map. For the first three orders of the MHWF, 
we have found amplifications greater than one for some scales. The same did not occur when we convolved the 
maps with gaussian filters taking into account the acceptance of the detector. In this case, the amplification 
for gaussian filters were less than one and in the best case, 10% greater than one. However, if we reject the 
acceptance of the detector, the gaussian filter provides amplifications greater than one for some angular scales 
in all the cases. In figure [5] (right) we expect that the maximum amplification would be obtained for a gaussian 
filter with dispersion equal to the one from the simulated source, since this is the scenario of a source embedded 
in a white noise background. This can be clearly seen in this figure. In figure [6] (right) we have a dipole with 
an excess of events 7% greater than the background in the direction of the galactic center, without acceptance. 
Due to the effect of the imposed dipolar anisotropy pattern, the gaussian filter does not provide the maximum 
amplification of the signal-to-noise ratio. In the case of the second dipole without acceptance, shown in the right 
side of figure [7J the gaussian filter provides the maximum amplification, since the simulated dipolar pattern is 
very faint (0.1%). In the case with sources and without acceptance, shown in figure [8] (right), the gaussian filter 
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Scale (degrees) Scale (degrees) 

Figure 6: Amplification for the dipole 1 case with (left) and without (right) acceptance. 




Scale (degrees) Scale (degrees) 



Figure 7: Amplification for the dipole 2 case with (left) and without (right) acceptance. 
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Figure 8: Amplification for the simulated sources case with (left) and without (right) acceptance. 
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Table 1: A7 (measured in degrees) for different anisotropies in the case of N src events from the source and an 
acceptance law for the detector. 
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Table 2: A7 (measures in degrees) for different anisotropies in the case of N src events from the source and no 
acceptance law. 
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provides an amplification greater than one, but still less than the ones provided by the MHWF. 

At first glance, the conclusion that the MHWF has a greater power of discrimination for sources when the 
background is not white noise seems reasonable. However, let us notice that the simulated source has a dispersion 
of 2.0°, and the peaks of amplifications are shifted to the right with respect to this value, in the case of the 
MHWF. So, the scale of the filter to be used and the scale of the source are not the same. Let 7 be the scale of 
the filter, corresponding to the dispersion of the gaussian of the used filter and 70 the dispersion of the simulated 
source, which in our case is 2°. We can define A7 = I7 — 7o| as being the distance between the scale of maximum 
amplification and the dispersion of the simulated source. This is shown in tables [l] and [2] . 



We notice in table [T] that A7 for different anisotropies is approximately constant for all filters, despite the 
anisotropy pattern of the background. However, for 50 events coming from the direction of the source, we notice 
that A7 fluctuates up to 1.6° in the case of the gaussian filter and at most 0.6° for the MHWF filters. The 
same is not true for the full sky simulations, presented in table [2j In this case, A7 fluctuates up to 1.4° for the 
MHWF filters and 0.8° for the gaussian kernel. 

6 Conclusions 

In this work we have considered a sinOcosO law for the zenith angle distribution for a virtual detector located 
at the same place as the Pierre Auger Observatory and an ideal detector observing the whole sky, with uniform 
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acceptance. We have simulated a source with different number of events coming from it, embedded in several 
backgrounds. We have calculated the amplifications of the signal-to-noise ratio in each case and verified that 
in the case of a source with dispersion 70 embedded in a white noise background, the maximum amplification 
is achieved with a gaussian filter with dispersion 7 = 70, as predicted. In the case of a faint dipole with no 
acceptance, the gaussian provided a greater amplification of the signal-to-noise ratio. However, if the background 
does not have a constant power spectrum, as in the case of the other anisotropy patterns, the Mexican Hat 
Wavelet Family filters provided, in general, a greater amplification of the signal-to-noise ratio, allowing one to 
identify point sources in maps containing arrival directions of cosmic rays. 

It is interesting to notice that the existence of an acceptance for the detector affects the power of discrim- 
ination of the gaussian filter. The amplification, in this case, is always below 1.1, whereas in the case of the 
MHWF the amplification of the signal-to-noise ratio can be more than ten. 

The acceptance of the experiment affects the power spectrum of the background in celestial maps. If it is not 
white noise, the maximum amplification of the signal-to-noise ratio is achieved by the MHWF kernel. However, 
even though we can amplify the signal-to-noise ratio using MHWF, we lose directional information about the 
associated cosmic ray events. If wc consider a gaussian filter and a white noise background, the optimal gaussian 
filter would have exactly the same dispersion as the source and the uncertainty on the position would arise from 
the resolution of the detector, rather than the analysis method used. By using MHWF, the greater the order of 
the wavelet, the greater the difference between the sizes of the source and the wavelet. Therefore, despite the 
gain on the power of discrimination of the source, there is a loss of resolution. 

Since the acceptance introduces a non-white noise background, for analysis involving the whole sky, such as 
blindsearcrj^J the MHWF amplifies the signal-to-noise ratio more than the gaussian. However, for small scale 
analysis, smoothing the maps with a gaussian kernel and reducing the area of scan could provide a greater 
amplification, specially if the background within the window of scan has a uniform power spectrum. 

A limitation to this technique is the projection of the celestial sphere to the plane. This approximation is 
good enough for small scales, when we can use the approximation sinO s» 6. For 9 > 7°, this approximation 
introduce bias and this method would have to be adapted to work on a spherical manifold. Wavelets on the 
sphere have been used on CMB studies [TP! ITT1 H"2l [T5] and presented good results. The next step of this work is 
a similar analysis using wavelets on the sphere, which allow us to search not only point sources, but also large 
scale structures in maps containing arrival directions of cosmic rays. 
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